Quantum Interference in Single Molecule Electronic Systems 



O 



R. E. Sparks\ V. M. Garcfa-Suarez^^^, D. Zs. Manrique\ and C. J. Lambert^ 

^ Department of Physics, Lancaster University, Lancaster, LAI 4YB, U. K. and 
Departamento de Fisica, Universidad de Oviedo & CINN, 33007 Oviedo, Spain 
(Dated: January 14, 2013) 

We present a general analytical formula and an ab initio study of quantum interference in multi-branch 
molecules. Ab initio calculations are used to investigate quantum interference in a benzene- 1,2-dithiolate (BDT) 
molecule sandwiched between gold electrodes and through oligoynes of various lengths. We show that when a 
point charge is located in the plane of a BDT molecule and its position varied, the electrical conductance exhibits 
a clear interference effect, whereas when the charge approaches a BDT molecule along a line normal to the plane 
of the molecule and passing through the centre of the phenyl ring, interference effects are negligible. In the case 
of olygoynes, quantum interference leads to the appearance of a critical energy Ec, at which the electron trans- 
mission coefficient T{E) of chains with even or odd numbers of atoms is independent of length. To illustrate the 
underlying physics, we derive a general analytical formula for electron transport through multi-branch structures 
and demonstrate the versatility of the formula by comparing it with the above ab-initio simulations. We also 
employ the analytical formula to investigate the current inside the molecule and demonstrate that large counter 
currents can occur within a ring-like molecule such as BDT, when the point charge is located in the plane of the 
molecule. The formula can be used to describe quantum interference and Fano resonances in structures with 
branches containing arbitrary elastic scattering regions connected to nodal sites. 
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INTRODUCTION 

The field of molecular electronics iH is a rapidly expanding 
research activity, which biidges the gap between physics and 
chemistry. Recently there has been much interest in develop- 
ing strategies to control the current through a single molecule 
lUlstl • Of the various effects that can be exploited, quantum 
interference is expected to play a fundamental role in long 
phase-coherent molecules |4], where multiple reflections can 
occur and in molecules made of rings, where electrons can 
follow multiple paths between the electrodes fslS- The mod- 
ification of the electronic properties of such systems has ap- 
plications such as the quantum inteiference effect transistor 
(QuIET) 101 and can potentially to be used for implementing 
data storage |8], information processing |j9|] and the develop- 
ment of molecular switches 1 10]. 

In this article, we study quantum interference effects in 
molecules between metallic leads using a combination of an 
analytical model and large-scale ab-initio simulations. We de- 
rive a versatile analytical formula for the electrical conduc- 
tance of molecular structures, which captures quantum inter- 
ference effects in linear and multi-branch molecules. For lin- 
ear oligoyne molecules or an atomic chain linking two elec- 
trodes, we predict that for odd or even-length chains, quan- 
tum interference leads to the presence of a critical energy Ec, 
at which the electron transmission coefficient becomes inde- 
pendent of length for odd or even numbers of atoms in the 
chain. The presence of this critical energy in more realistic 
structures is confirmed by performing an ab initio calcula- 
tion of electron transmission through an oligoyne molecular 
wire connected to gold electrodes. We also present results of 
an ab-initio numerical simulation on an electrostatically-gated 
benzene dithiol (BDT) molecule, attached to gold electrodes, 
which is an example of a QuIET. In this calculation, gating is 



achieved through the presence of a calcium or potassium ion, 
which induces quantum interference as the position of the ion 
and the molecular orientation are varied. We show that the 
qualitative features of this interference effect are captured by 
the above analytic formula through an appropriate choice of 
parameters. Finally, we note that quantum interference in such 
multi -branch structures leads to the appearance of large inter- 
nal counter currents, which exceed the external current carried 
by the electrodes. 



AN ANALYTICAL FORMULA FOR ELECTRON 
TRANSPORT THROUGH MULTI-BRANCH STRUCTURES 

Fig. ([T]) represents a tight-binding (Hiickel-type) model of 
a multi-branch structure, in which each atom is assigned a sin- 
gle atomic orbital. The structure consists of left and right leads 
connected to external electron reservoirs (not shown). The 
atoms of the left lead (L) are labelled j = —1, —2, —3, .... 
The orbital energy of each atom is denoted el and these 
are coupled to each other by a nearest-neighbour matrix el- 
ement — 7l. Similarly, the atoms of the right lead labelled 
j = 1, 2, 3, ... , are assigned orbital energies and these 
are coupled to each other by a nearest-neighbour matrix el- 
ement —"fR. The loop structure comprises AI branches, la- 
belled I — 1,2, ... , M. Branch / possesses Ni atoms, labelled 
ni = 1,2, Ni, with atomic-orbital energies coupled by 
nearest neighbour matrix elements —7/. (Note that hopping 
matrix elements could be positive or negative and the inclu- 
sion of a minus sign is merely convention. For simplicity, we 
consider the case of a real hamiltonian, since in molecules, 
orbital effects due to applied magnetic fields are usually neg- 
ligible.) The left-most atom (n/ = 1) of each branch is con- 
nected by a matrix element —a/ to a nodal atom (labelled L) 
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FIG. 1: A multi-branch structure with nodal sites L and R on the left and right, connected to external current-carrying leads, by hopping matrix 
elements — ql on the left and —Pr on the right and to internals branches {?} by hopping matrix elements —a; and — /3; respectively. The 
energies of the nodal sites are e'^ and e'^ . The site energy and hopping matrix element of branch I are e; and —7; respectively. 



of orbital energy e° . The latter is connected to the right-most 
atom of the left lead by a matrix element — a^. Similarly, the 
right-most atom (n; = Ni)of each branch is connected by a 
matrix element — to a nodal atom (labelled R) of orbital en- 
ergy e^f, which in turn is connected to the left-most atom of 
the right lead by a matrix element 

In the presence of an incoming plane wave from the left, 
the solution to the Schrodinger equation, H-^} — Eip, in the 
left lead (j < —1) is of the form 

Similarly, the solution in branch I can be written 

= Aie''"'" + Bie-''"'" (2) 
and the wavefunction in the right lead {j > 1) is of the form 

^pf) = t{E)e'^'-'^ (3) 

Finally, the wavefunction on the left and right nodal atoms 
will be denoted xl and xr respectively. In the above equa- 
tions, E is the energy of the incident electron and r{E) and 



t{E) are transmission and reflection amplitudes. For a given 
E, the dimensionless wavenumbers in the left and right leads, 
and in branch I are given by fc,, = cos^^ (^^-T^) ^^ere 
the index 77 is either L, R 01 1 respectively. The correspond- 
ing group velocities {a,f/h)dE/dk can be written (a^/?i)u^, 
where a^j is the atomic spacing in region 77, w,, — 27,, sinfc^. 
In what follows, we adopt the convention of choosing real val- 
ues of kjj, such that is positive and complex values of fc^, 
such that Im(A;,,) is positive. 

Our initial goal is to obtain an expression for the transmis- 
sion amplitude t{E), which as shown in the appendix, can 
be obtained either by matching wavefunctions at the nodal 
atoms or by using Green's functions. According to the Lan- 
dauer formula, the zero-bias electrical conductance is simply 
{2e'^/h)T{Ep), where Ep is the Fermi energy and 

T{E) - {vr/vlME)^, (4) 

which satisfies T{E) + R{E) = 1, where R{E) = \r{E)\'^ 
is the reflection coefficient. In terms of T{E), the current per 
unit energy carried by the left and right leads is {2e / h)T (E) 



3 



and since T{E) < 1, the current per unit energy in the left 
and right leads cannot exceed 2e/h. As we shall see below, 
for M > 1, this upper bound does not apply to the current per 
unit energy carried by the internal branches, which we denote 
{2e/h)Ii. Indeed for M > 1, // can be either positive or 
negative and is unbounded. 

As shown in the appendix, T{E) can be written 



T{E) = VL 



\G 



RL\ 



IB. 



VR 



(5) 



This expression is very general and shows how the various 
contributions combine to control the current through a single 
molecule. Equation (|5j shows that the transmission coeffi- 
cient T{E) is a product of several factors; the "group veloci- 
ties" vl and describe the ability of the left and right leads 

to carry a current, (^^^ and (^^^ describe the ability of 
the couplings between the nodal atoms and the external leads 
to transfer electrons and finally Grl describes the ability of 
a current from a source at node L to be carried to a current 
sink at node R. In this expression, Gbl describes propagation 
from the nodal site L to at the nodal site R and is sensitive to 
quantum interference within the multi-branch structure. Since 
Vl and vji have dimensions of energy, whereas Gbl has di- 
mensions of energy the right hand side of equation (|5]l is 
dimensionless, as expected. 

As shown in the appendix, Grl is given by 



'-'RL — 



where 



A = - (oL - XL){aR - xr). 
In this equation. 



xf = a^i siiifc;(7Vi)/[7isinfci(7V/ + l)] 



and 



xf = sinfci(^0/[7i sinfcK^i + 1)] 
Finally, the parameters and are given by 



L Akr 



aL = {el -E)-^e 

IL 



(12) 



(13) 



(14) 



and 



aR 



(4 



E) - ^e*''^ 
IR 



(15) 



Clearly the parameters and aR are independent of the de- 
tails of the internal branches I and are properties of the left and 
right leads and their respective nodal atoms only. Properties 
of the branches are contained within the parameters xl, xr 
and y only. From Eq.(|6]l, T{E) will vanish when y=0. This 
condition for destructive interference does not depend on the 
parameters describing the leads (e\, e\, 7l, 7fl). Nor does it 
depend on the parameters describing the contacts to the leads 
(oLh, £L,otR, £r)- It is a fundamental property of the branches 
and their couplings to the nodal sites. 

As noted in the appendix, equation (|5]) is extremely gen- 
eral. With a slight modification of the nodal energies e° and 
e^j, it can be used to describe the effect of Fano resonances 
(6) due to dangling bonds at the nodes. Furthermore, with a slight 
redefinition of yi, xf and xf", it describes electron transmis- 
sion arising when the branches are replaced by arbitrary elas- 
tic scatterers connected by single bonds to the nodal sites. 
An alternative form of equation (|5) is obtained by writing 
^"^-^ A = Ai + iA2, ai, = gl — if l and aR — clr — iTr, where 
CLL = el - E - {al/jL) cos fc^ and Tl = (ai/7i) sin fci 
and similarly for clr and Tr. With this notation. 



M 

y^^yi, 
1=1 



(8) 



Ai = y^ - {xL - aL)ixR - gr) + rLtR (16) 



and 



and 



where 



M 

XL=^xf 
1=1 



M 

XR^^X^ 

1=1 



(9) 



(10) 



yi = aiPi sinfci/[7i smki{Ni + 1)] 



(11) 



A2 ^Tr{xl - gl) +TLixR - gr) (17) 



and 



T{E) = 



^TLTRy^ 



A? 



A2 
^2 



(18) 



Equation (|5]l describes the transmission coefficient of the 
combined structure and allows us to evaluate the current per 
unit energy {2e/h)T{E) due to incident electrons from the 
left lead with energies E. We shall also be interested in the 
current per unit energy (2e//i)/; carried by branch I. As 
shown in the appendix, this is given by 
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FIG. 2: A single impurity in a 1-d ciiain 



FIG. 3: A schematic of the nodal atoms connected by an atomic 
chain. 



/, - T{E)yi/y, 



which clearly satisfies 



M 



(19) 



(20) 



1=1 



Unlike T{E), which satisfies < T{E) < 1, Ii can have 
arbitrary sign and arbitrary magnitude. 

Before using equation (|5]l to describe quantum interference 
within linear and multi-branch molecules, we consider the 
simplest choice of a single impurity level, weakly coupled to 
external left and right leads, by matrix elements ai and /3i 
respectively, is shown in Fig®. 

This corresponds to the choice M = 1, Ni = 1, jl = 
IR = l^R, £l = £r = £% In this case equation (|5]l 
reduces to the well-known Breit-Wigner formula 



T{E) = 



(21) 



where ctl = (q!?/7l) cos fc^, aR = {/Sy^n) cos kR, Ti 
("i/7l) sinfci, andTij = {l3f /jr) sin kR. 



QUANTUM INTERFERENCE IN LINEAR MOLECULES OR 
ATOMIC CHAINS. 

The choice Af = 1 A^i > 1 corresponds to the case of exter- 
nal left and right leads, coupled by matrix elements a l and 
/3r respectively, to nodal sites L and R, which in turn are con- 
nected by matrix elements ai and /3i to an atomic bridge of 
A'^i atoms. The case A^i = 3 is shown in Fig©. 
For AI — 1, one obtains 



sin ki 



7i sinfci(A^i + 1) 



(22) 



XL 



xr 



sin fciTVi 



71 sinfci(A^i + 1) 



Pi sinfciTVi 
71 sinfci(7Vi + !)■ 



(23) 



(24) 



In the case of a metallic or "tt bridge", fci will be real. In 
the case of a "a bridge", (which acts as a tunnel barrier), ki 
will be imaginary and equation (|5]l (or equivalently equation 
(fTsT l) describes electron transport via superexchange. Equa- 
tions (l22l i. (l23T l. (l24l i highlight a curious feature, which oc- 
curs at a special energy Ec, which corresponds to electrons 
propagating at the band centre of a vr bridge and at which 
fci — 7r/2. At this energy, y, xl and xr become indepen- 
dent of the length A^i of the bridge. On the one hand, if the 
bridge contains an even number of atoms (ie if A^i is even), 
then XR^XL^Q,y = (ai/JiM)!-!)'"^^/'^ and 



AVLTRiaipihif 



{{aiPihi? - aL~aR + TlTr)^ + {dLTR + aRTL^ 

(25) 

On the other hand, if the bridge contains an odd number of 
atoms, then xl, xr and y diverge and 



TiE, 



(26) 



which is independent of the length A^i of the bridge. This 
situation can arise, for example, in the case of oligoynes con- 
nected to external electrodes. 

These predictions are shown in Fig.(|4) for increasing num- 
bers of atoms in the wire, Ni = 2, 4, 6 and 8. At the criti- 
cal energy Ec ~ 0.5eV, all curves intersect. Consequently, 
for energies E slightly greater than Ec, T{E) will either in- 
creases monotonically as the length of the bridge increases by 
2, and for E slightly less than Ec, T{E) will decrease when 



-2 2 

E-Ep[eV] 

FIG. 4: Transmission functions for increasing lengths of molecular 
wire using the general formula. The system is modelled using the 
following parameters, in the leads el = er = 0.0 and = ~ 
5.0, for the contacts El = = —0.4, in the chain ei — 0.5 and 
7i — 6.0 and the coupling between the wire and the electrodes are 
defined by q:l = Pr = 3.0. All curves intersect at Ec ~ 0.5eV. 
Close inspection reveals that E ~ — 0.2eV the curves approach each 
other, but do not intersect at a single energy. 
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FIG. 5: Transmission curves from the SMEAGOL simulations of 
oligoynes of varying lengths 

the length of the bridge increases. This effect is a clear mani- 
festation of phase-coherent quantum transport. 

To demonstrate that this effect is present in atomistic calcu- 
lations of electron transport, we compare equation (l25i with a 
calculation based on the ab-initio transport code SMEAGOL. 
This code uses a combination of density functional theory 
(DFT) 1 12] and the non-equilibrium Green's function formal- 
ism 1113(1 to calculate the transport characteristics of atomic 
scale devices. The DFT Hamiltonian is obtained from the 
SIESTA code Q and is used by SMEAGOL to calculate the 
electronic density and the transmission. Within the NEGF the 
system is divided in three parts, the left lead, the right lead and 
the extended molecule (EM). The EM contains the molecule 
plus some layers of gold, whose electronic structure is mod- 
ified due to the presence of the molecule and the surfaces 



(a) (b) 

FIG. 6: Molecular structure used in the transport simulations with 
the potassium point charge in configuration CI (a) and the calcium 
point charge in configuration C2 (b). 

and differs from the bulk electronic structure. The molec- 
ular structure consists of an oligoynes capped with phenyl 
rings and attached to the electrodes by thiolate groups. The 
SMEAGOL results are shown in Fig.©, which clearly pos- 
sesses a critical energy Ec ~ 0.5eV at which all curves (at 
least for the longer chains) intersect. The analytic expression 
assumes that the parameters ei and 7 describing the chains are 
independent of length. In fact the self-consistent DFT param- 
eters of the shortest chain (A^i = 2) differ slightly from those 
of the longer chains and therefore the black curve of figure 5 
does not quite pass though the intersection point at E = Ec- 

Clearly the length independence of even and odd chains 
leads to an even-odd oscillation in the electrical conductance 
of oligoynes, when Ef is close to Ec- This effect has also 
been observed in experiments on atomic wires of Au, Pt, and 
Ir IitIi . which exhibit electrical conductance oscillations as a 
function of the wire length and similar oscillations as a func- 
tion of bias voltage and electrode separation [[IslH. Several 
theoretical papers ll20ll -||34ll have also addressed these osilla- 
tions. The above analysis also demonstrates that this effect is 
present in multi-branch structures, provided the band centres 
of different branches occur at the same energy. 

QUANTUM INTERFERENCE IN A TWO-BRANCH 
MOLECULE 

We now turn to the quantum interference effect transistor 
(QuIET) discussed in |7], which corresponds to the choice 
M = 2. To demonstrate that equation (|5]l (or equivalently 
(fTSb ) reproduces the key features of a QuIET, we compare 
it with the results of a detailed simulation using SMEAGOL 
Ijjfl. 

The atomic arrangements for the SIESTA/SMEAGOL cal- 
culations are shown in Fig. (|6]l. The first arrangement (CI) 
corresponds to the point charge located along a line perpen- 
dicular to the plane of the molecule and which passes through 
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FIG. 7: Transmission functions obtained from SMEAGOL simula- 
tions with the potassium point charge moving closer to the molecule 
through (a)-(d). The continuous curve corresponds to CI and the 
dashed curve to C2. 



FIG. 8: Transmission functions obtained from SMEAGOL simula- 
tions with the calcium point charge moving closer to the molecule 
through (a)-(d). The continuous curve corresponds to CI and the 
dashed curve to C2. 



its center. In this configuration the point charge produces a 
symmetric voltage which affects the two branches to the same 
extent. The second arrangement (C2) corresponds to the point 
charge located in the plane of the molecule, closer to one 
branch of the BDT. In this case the two branches are subject to 
different electrostatic potentials, which induces quantum in- 
terferences in the electron transmission through the molecule. 
Both configurations were simulated using a point charge of 
either potassium (K) or calcium (Ca), giving a total of four 
cases. K and Ca are alkali and alkaline-earth atoms with 1 
and 2 valence electrons in the last shell, respectively. Due 
to their high electropositivity both atoms lose their valence 
electrons when they are inserted in the unit cell and become 
ionized with a charge of +e and H-2e, respectively. The com- 
plete removal of the valence electrons from these atoms can 
be ensured by reducing the cutoff radii of their orbitals to 3.5 
Bohr, which confine the electrons in the atom more closely 
and therefore increase their energy, making sure they move 
to lower energy states in the extended molecule. The ba- 
sis sets used in the simulation were single-zeta (SZ) for the 
point charge and double-zeta polarized (DZP) for all other el- 
ements. The exchange and correlation potential was calcu- 
lated with the generalized gradient approximation (GGA) and 
the Perdew-Burke-Ernzerhof parametrization [15|. The gold 
leads were grown along the (001) direction, and each side of 
the extended molecule had 3 and 5 layers, respectively, with 
36 atoms (12x3 atoms) per layer. The molecule was con- 
tacted in a hollow configuration to four additional gold atoms 
on each side. Since the system was much larger in the x (3 
atoms) than in the y direction (12 atoms to leave space for the 
charge to move), 1 fc-point was used along x and 4 A: -points 
along y. 

The results are shown in Fig. d?) and Fig. (|8]l, for potassium 
and calcium, respectively. Each graph contains two curves, 
corresponding to the cases CI and C2. In plot (a) the charge 
is located at a far distance, 15 A from the molecule, and 



therefore both CI and C2 produce the same curve. From (b) to 
(d) the charge is gradually moved towards the molecule (6.29, 
5.29 and 4.79 A away from the center of the ring in (b), (c) 
and (d), respectively). 

We observe that when the charge moves towards the 
molecule the peaks shift in energy in the negative direction 
due to the positive potential. However the effect is different 
depending on where the charge is located relative to the ring. 
As can be seen, there is a clear difference in both Fig. ([7} 
and Fig. dHJ between the continuous and dashed transmission 
curves in graphs (b)-(d). An extra peak in the dashed trans- 
mission curve (C2) appears and the height of the HOMO peak 
is reduced, whereas the continuous transmission curve (CI) 
is simply shifted to lower energies without much change in 
the resonances. Also, through comparison of (a) and (d) we 
notice a clear narrowing of the HOMO and broadening of the 
LUMO peak in all cases. We observe a clear reduction of the 
transmission at the Fermi energy when the charge is located 
closer to one arm of the molecule (C2). In contrast, for sys- 
tem CI, there is very little change of the transmission about 
the Fermi energy, because the point charge produces the same 
phase shifts in the two branches and therefore does not modify 
interference effects associated with coherent superposition of 
waves propagating along separate paths. 

We also checked the projected density of states (PDOS) on 
each branch of the BDT to see the specific effect of the charge 
on the electronic structure in each case. In CI the PDOS on 
each branch remains equally distributed and simply shifts to 
lower energies. However, in C2 there is a clear difference in 
the PDOS on each branch; the PDOS on the closest branch to 
the charge is more affected and shifted to lower energies than 
the PDOS of the opposite branch. This supports the observa- 
tion of the previously suggested QIE. 

To elucidate the underlying physics, we employ Eq. (|5]l to 
model electron transmission through a two-branch structure. 
In the absence any charge, we choose the hopping parameters 
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FIG. 9: Transmission curves for the tight binding model. The contin- 
uous hne corresponds to the case where the charge is located along 
a line perpendicular to the ring and which passes through its center. 
The dashed line corresponds to the case where the charge is closer to 
one arm. Panel (a) shows the transmission curves when the charge 
is away from the molecule. From panel (b) to (d) the onsite energy 
parameters are changed to simulate the charge moving towards the 
molecule. 



= i^R — 1-5, 7l = IB. = 4.0 and 71 = 72 = 1.0 and 
the on-site energies £^ — £% — 2.0, ei = 0.0, £2 = 0.0 and 
El — £r — 4.0. This leads to the transmission curve shown in 
Fig. ^ (a), which is very close to the ab-initio result. In con- 
figuration CI, where a charge affects both branches equally, 
the presence of a charge is modelled by shifting the on-site en- 
ergies El, £1, £2 and er, upwards or downwards by the same 
amount, which depends on the sign and strength of the charge. 
The outcome produced by a positive charge is represented by 
the continuous transmission curves in Fig. (|9]i. The charge 
moves closer to the ring from (b) to (d) and the parameters 
are chosen as follows (b) £° = £5; = 1.4, £1 = £2 = —0.8, 



(c) £' 



_ 

L 



= 1.2, £1 = £2 = -1.8, (d) eI= E% ^ 1.0, 



£1 = £2 = —2.0. In each of these plots, El — er remain 
unchanged throughout. As in the ab-initio simulations we see 
that the entire transmission curve is shifted to lower energies 
and quantum interference effects are negligible. Interestingly, 
as a consequence of this shift and the corresponding change 
in the electronic structure, the width of the variability in the 
local density of states at the contact, the width of the HOMO 
decreases and the width of the LUMO increases, in agreement 
with the ab-initio results. 

To produce quantum interference, we now examine the ef- 
fect of a scanning point charge placed in configuration C2; 
i.e. closer to one branch of the ring. To model this effect us- 
ing Eq. (fTSl l. the parameters are now adjusted asymmetrically; 
i.e. they are changed less in the branch which is far away from 
the charge and more in the branch which is closer. The adjust- 
ment also includes changing the contact points £^ and e'j^ as 
these will feel a smaller effect from the charge than the nearer 
branch. The adjusted parameters are chosen as follows (b) 
£l=£0^ = 1.35, £1 = -1.0, £2 = -0.5 (c) eI=e°j^^ 1.1, 
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FIG. 10: (a) Total current through the molecule, (b), and the current 
through branch 1, (c), and branch 2 of the molecule. Parameters used 
are exactly as the dashed red line in Fig 9(d), namely er = el = 
4.0, el = e'A = 1-0, ei = -1-8, £2 = -0.9 7^=7^ = 4.0, 
71 = 72 = 1-0, aL = I3r = 1.5. 



= -1.5, £2 = -0.8 (d) £« = E% = 1.0, £1 = -1.8 
and £2 = —0.9. As before, el = er and are unchanged. 
The transmission corresponding to these parameters is shown 
Fig-® (dashed curves), where the point charge is brought suc- 
cessively closer to the molecule from (b) to (d). We see again 
from (a) through to (d) that the peaks have all shifted to lower 
energies, but the HOMO dramatically changes and reduces 
its height. Also, an additional peak appears due to the point 
charge effect on the electronic structure on only one arm of 
the molecule, which causes interferences in the transmission 
through the system. This again agrees with the SMEAGOL 
simulations and suggests this analytical model captures the 
qualitative features of transmission in ring-like molecules. 

Having established that the analytical model captures the 
essential features of the ab-initio simulations, we now show 
how this model can be employed to examine the internal cur- 
rents within different branches of the molecule, which are ob- 
tained by evaluating equation iT% . When the ion located close 
to branch 1, the lower graphs of Fig. (fTOl l show the internal 
currents Ii and I2 through the individual branches, whereas 
the upper graph shows the total current /tot = Ii + I2 = 
T{E). Figs. ( [Tot (b) and (c) clearly demonstrate that the cur- 
rent in a single branch can greatly exceed the total current 
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through the molecule when a counter current of opposite sign 
occurs in the other branch of the molecule and can clearly ex- 
ceed the upper bound of /tot ~ 1- The appearance of such 
unbounded counter currents is yet another manifestation of 
quantum interference within single molecules. I,16J . 



SUMMARY 

In conclusion, we have presented ab initio simulations and 
an analytical formula, which highlights a range of interfer- 
ence effects in single and multi-branch structures. The ana- 
lytical solution is rather versatile and has the advantage that 
it can be evaluated on a pocket calculator. It provides insight 
into length-independent electrical conductances for even and 
odd oligoyne chains, when the Fermi energy coincides with 
the band centre of the oligoyne bridge and allows us to pre- 
dict that this behaviour is also present in multi-branch struc- 
tures, provided the branches share a common band centre. As 
demonstrated in the manuscript the energy Ec at which this 
odd-even effect occurs corresponds to kl = Tr/2. This con- 
dition is very general and is independent of the nature of the 
orbitals. For the particular case of oligoynes, this is a conse- 
quence of TT coherent transport, but for other systems, such a 
metallic wires, this would not be the case. The case M = 1 
demonstrates that quantum interference does not require the 
presence of physically different paths, because even in this 
case, interference due to scattering from nodal impurity sites 
and connections to external leads arises from the amplitudes 
aL and qr in equations (14) and (15). Both the magnitudes 
and phases of these amplitudes appear on the right hand side 
of equation (7) and therefore even for a single-branch system, 
they contribute to interference. 

Ab initio simulations based on density functional theory, 
demonstrate the presence of quantum interference in BDT, 
due to electrostatic interactions associated with a scanning 
point charge positioned close to the molecule. We have shown 
that a scanning charge located within the plane of a BDT 
molecule produces a sizeable quantum interference, whereas 
a charge approaching the molecule along a line perpendicular 
to the plane produces a much smaller effect, in agreement with 
the analytical formula. In spite of the consistency between the 
TB result and the ab initio result for the BDT system, there are 
of course quantitative differences between them. In part this 
arises because the tight-binding model includes only a sin- 
gle ("tt") orbital per atom, whereas the ab initio description 
includes both tt transport and a tunneling. In addition, the 
tight-binding model includes only a single scattering channel 
in each lead, whereas the ab initio model contains multiple 
channels. 

Using the analytical model, we have also investigated the 
internal currents within a two-branch molecule and demon- 
strated that large currents and counter currents can occur, 
which exceed an upper bound for the total current through the 
molecule. 

We wish to thank the Spanish Ministerio de Ciencia e Inno- 



vacion, the UK EPSRC and the European Research Networks 
NanoCTM and FUNMOLS for funding. 



APPENDIX I: DERIVATION OE EQUATIONHFOR 
TRANSMISSION THOUGH THE MULTI-BRANCH 
STRUCTURE OF FIGURE (1). 

We derive the equation for T{E) by matching wave func- 
tions at the nodes of a multi-branch structures and later make 
a comparison with results obtained from a corresponding 
Green's function analysis. The starting point is the tight bind- 
ing Schrodinger equation, which can be written 



(27) 



where the summation is over all nearest neighbours v of site 
II. Choosing /i to label the site just to the left of the nodal 
site L (whose wave function is denoted xl) yields e^V^i^i' — 

lL->p''-2 - c^lXl = Eil:^^}, where 7^ = 7-I, -2 and ul = 
7l _i. From this expression, and noting that the Schrodinger 
equation in the left lead takes the form of a recurrence relation 
[i35il . the wave function at the node L is given by 



XL 



lL{l + r) 



CtL 



(28) 



Similarly, choosing to label the site just to the right of the 
nodal site R (whose wave function is denoted xr) yields 



XR 



iRt 

Pr ■ 



(29) 



Choosing fj, to label the first site (1;) of chain I yields for all I, 

XL = {li/ai){Ai+Bi) (30) 

and choosing fi to label the last site {Ni of chain I yields for 
all I, 

XR = (7;//3/)(^;e*W+i) + Bze-^W+i)) (31) 
Finally choosing ji to label the nodal sites L and R yields 

M 

4XL - "'^i'^ - + KS)e''" ) = EXL (32) 



1=1 



and 



M 



SrXR - ^'^S - f^RtiEy"" - ExR (33) 



1=1 



Equations ( l28T l with ( |32] | and(|29l) with ( l33T l combine to yield 
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XL-^(l + r) 



(0 



2aii sin(fci) 



and 



(34) 



aR 



(35) 



where and an are given by equations (fl4l i and (flST l. 

From the form of the wave functions in the branches, given 
by equation (2), these can be written 



XL 



1 / —2iaL sin(fci) 
^ I 



(36) 



M 

E 

1=1 



aL 




Since equations dSOl l and dJTT i yield 



This is demonstrated by noting that the Green's function 
for a finite Unear chain of Ni sites, with nearest-neighbour 
hopping elements —7/ and diagonal elements ei is 



. ,s _ Ai sin kim sin ki[n'i - (Ni + l)] (fern/ < nj) 
9i{ni,ni)- ^jgin/fc,n;sin/c/([n/ - (TV, + 1)] (forn;<n/) 

(41) 

where Ai = l/[7; sin A:/ sin ki{Ni + 1)]. An alternative form 
of this expression is 3/ (n;, nj) = {Ai / 2) {cos ki[Ni + l — \ni — 
n[\] -coski[Ni + l-ni-n[]) 

The quantity gi{ni,n[) is the Greens function matrix ele- 
ment connecting atom n; to atom nj of the decoupled branch 
I, which would arise when ai — f3i — 0. The off-diagonal 
matrix element describing propagation from one end of such 
a branch to the other is 



5/(1,, TV,) = -sinA://7,sinA:/(iV, + 1) 



(42) 



whereas the diagonal matrix element evaluated on an end atom 
is 



A, 



1 



aie 



Hki(Ni + l) 



XL 



Bi ) -2^7^sinfc^(iV^ + l) V-aie*'='(^'+i) A ) \xR. 

(37) 

Ai and Bi can be eUminated from equation (|36] | to yield 



5/(1,, 1,) = 5,(iV,,iV,) = -sinA:/iV,/7,sinA:/(iV, + 1) 

(43) 

s expected, these quantities diverge when sin ki (7V; + 1) — 0, 
which corresponds to the eigenenergies of an isolated branch. 
In terms of these Greens functions. 



XL 
XR 



G 



2iaL sin(fci) 




LL 
RL 



In this expression, the matrix G has the form 



G = 



Gll Glr 

Grl Grr 



2iaL sin{kL) 
(38) 



(39) and 



yi = -ai(5igi{l,Ni) 



(44) 



(45) 



and is given by 



xf = -P^gi[NuNi) 



(46) 



-aL 
-aR 



XL y 
y XR 



(40) 



where xl, xr and y are given by equations (|9]l and ( fTOl l. 
From this expression, one obtains XR and hence the transmis- 
sion amplitude t, via equation ( |29] |. 

The physical meaning of the various contributions to the 
above expressions can be understood by carrying out a par- 
allel analysis based on Green's functions |3^ 37J, which re- 
veals that equation ( l40b is simply Dyson's equation for the 
Green's function matrix elements involving the nodal sites L 
and R. Comparison with Refs. ll36[ 37 1 also demonstrates that 



VLiotL/lLY and VR{aR/^Rf in Eq,(5) are imaginary parts 
of the self-energies of the left- and right-hand electrodes, re- 
spectively. 



Within a Green's function approach, one defines the nodal 
self energy matrix a to be 



M 

fT = ^fT/ 

1=1 



(47) 



where cri is the contribution to the self-energy from branch /, 
given by 



-ai 
-/?/ 



91 



-Oil 
-/?/ 



(48) 



In this expression 5/ is the Greens function connected the end 
atoms of an isolated branch: 
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9i{h,h) 
9i{Ni,h) 



9iih,Ni) 
gi{Ni,Ni) 



(49) 



This demonstrates that 



XL y 
y XR 



(50) 



and therefore equation (l40l i takes the form of Dyson's equa- 
tion: 



G- 



9l 




9n' 



(51) 



where 



and gji — —a„^ are diagonal elements 



of the Greens function of the decoupled semi-infinite chains 
(obtained by setting all 0), evaluated on the left 

(L) and right (R) nodal sites respectively. This also demon- 
strates that the form of equation(|5) and in particular Gbl 
does not change even when the branches I are replaced by 
arbitrary elastic scattering regions, connected to nodal sites 
by bonds ai and Pi, provided gi is replaced by the Green's 
function of the scattering region. With this redefinition 
of yi, the condition for destructive interference {y = 0) re- 
mains unchanged. For example, if instead of a linear chain 
of ni sites, branch I is replaced by a loop of ni sites, , then 
equation ( 1411 1 is replaced by the Greens function of a linear 
chain of ni sites with periodic boundary conditions, namely 
gi{ni,n[) = {coski[Ni/2-\ni-n'i\])/{2-fsmkism[kiNi/2]) 
and equation ( |49] | is replaced by 
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gi{ni,ni) gi{ni,mi) 
gi{mi,ni) gi{mi,mi) 



(52) 



where ni and mi label the sites of the loop connected to the 
nodal sites L and R respectively. Taking this to an extreme, 
any of the branches I could even be replaced by a multi -branch 
scatterer, simply by replacing gi by the Greens function of 
an isolated multi-branch system, obtained from G by setting 
aL^I^R^ 0. 

The above analysis, which focusses on the wave-like nature 
of Greens functions is rather different in spirit from alternative 
approaches which emphasise the algebraic nature of Greens 
functions, which for finite structures, take the form of ratios of 
polynomials, whose denominator is proportional to the secular 
equation |38]. To make contact with this approach, we note 
that equation (ISTT i yields 



X - an -y 
-y X -UL 



(53) 



where A = Ai + zA2, and therefore the equation Ai = 
is the secular equation for the isolated multi-branch structure. 



which arises when aL = f^R = 0. More generally, from 
equations (fTST i and (fTTT i. the equation Ai — TlTr is the sec- 
ular equation for the same isolated system, but with the site 
energies of the nodal atoms shifted by the real part of their 
respective self energies. 

Finally the current per unit energy in branch I, carried by 
electrons of energy E injected from the left lead is (2e//i)/;, 
where 



Ii = ^{\Aif-\Bif) 



(54) 



Expressions for Ai and Bi are obtained from equation (l37T i. 
which combine to yield equation(fT9]l of the main text. 

The above comparison between the wave-function- 
matching and Green's function underpins a deep understand- 
ing of equation (1381 ). because if /i < — 1 labels a site in the left 
lead and v > n labels a site inside the scattering region or in 



the right lead, then the wave function V'l 
the expression 



is related to G^.^j. by 



(55) 



Furthermore, starting from the limit — Q and then in- 
cluding the effect of via Dyson's equation yields Gll — 
{lLlaL)e-'''^GL^-i and Grl = (7L/aL)e-''^'^Gfl,-i. 
Hence equation ( |38]) can be written in the intuitive form 



XL 
XR 



Gl-1 
Gr-1 



2z7Lsin(fcL), (56) 



which is simply an example of equation (1551 1. with /i = — 1 
and v = L 01 R. 

As mentioned in the main text, equation(|5]l is extremely 
versatile. For example, the case of M = 1 A^i > 1, can be 
used to describe a donor-bridge-acceptor molecules. In this 
case, to obtain a simple description of rectification, all param- 
eters should be assigned and appropriate dependence on the 
applied voltage V. The simplest model is obtained by setting 

;0(G) + el//2, e«(l/) = 



eV/2, and then computing 

EF+eV/2 , 



tL{V) = 6^(0) + eV/2, tl{V) : 
6«(0) - eV/2 e%{V) = e%{0) - 
the current via the expression / — j^^Z'^y'/2 T{E)dE 

To further demonstrate the versatility of equation([T9), we 
end this appendix by noting that it readily describes the effect 
of Fano resonances on transport and the effect of coupling to 
a molecule at different points along its length. To illustrate 
this, consider a structure in which dangling branches, labelled 
I — and I ~ M + 1, are attached by couplings and Pm+i 
to the nodal sites on the left and right respectively, as shown 
in figure (fTTT l. 

In this case, equations (|5]) and ( fTSl ) are unchanged, except 
that e^? and are renormalised by the self energies of the 
dangling branches, and replaced by 



-0 



£fl + P\i+i9m+i 



(57) 
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FIG. 11: A diagram to show the case of a molecule with dangling 
branches. The top and bottom panels show two equivalent represen- 
tations of the same system. 



and 



where 



"0.90, 



50 



sin fco^o 



7osinfco(^o + 1) 



(58) 



(59) 



and 



sin k 



9M+1 — 



m+iNm+1 



jM+i sin kM+iiNM+i 



1 



(60) 



Clearly an anti -resonance occurs when the energy E coincides 
with an eigen-energy of either of the two branches, because 
at these energies, one of the Green's functions go or gM+i 
diverges and therefore one of the renormalised nodal site en- 
ergies e^j or £^ diverges. This is equivalent to introducing 
an infinite potential at one of the nodes and therefore at these 
energies, T{E) vanishes. This behaviour arises from the in- 
teraction between bound states in the dangling branches and 



the continuum of states associated with the external leads and 
is typical of a Fano resonance. 

By redrawing Figure figure (ITTk ) as shown in figure 
(fTTb). one can see that the above equation describes a linear 
molecule contacted at atoms within the length of the molecule, 
rather than simply at the end atoms. As an example, con- 
sider the case when M — 1, Si — Eq ~ 62 ~ £l — £r, 
71 = 7o = 72 = cko = /?2 = "1 = The system then com- 
prise a linear chain of length L = A^o + -^1 + ^2 + 2 sites, 
connect to external leads by nodal sites located at positions 
A^o + 1 and A^i + A^o + 2 along the chain. By varying A^o ™d 
A^i but with fixed L, the expression for T{E) then describes 
quantum interference effects which arise when external leads 
are connected to a fixed length molecule, at different locations 
along its length. 
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